Rectification of light refraction in curved waveguide arrays 
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An 'optical ratchet' for discretized light in photonic lattices, which enables to observe rectification of light 
refraction at any input beam conditions, is theoretically presented, and a possible experimental implementation 
based on periodically-curved zigzag waveguide arrays is proposed. © 2010 Optical Society of America 
OCIS codes: 130.2790, 230.3120, 230.7370, 000.1600 



Transport of discretized light in photonic lattices has 
attracted a great attention in recent years from both 
fundamental and applied viewpoints [1,2]. Engineered 
photonic lattices offer the possibility of tailoring diffrac- 
tion [3] , refraction [4] , and of shaping polychromatic light 
beams [5-7] . They also provide a rich laboratory system 
to mimic the dynamics of quantum particles driven or 
disordered periodic lattices, with the observation of the 
optical analogues of Bloch oscillations [1,8], Anderson 
localization [9, 10] and dynamic localization [11-13] to 
name a few. Recently, the transport of ac-driven quan- 
tum particles in periodic lattices has received a great 
interest also for the possibility to achieve Hamiltonian 
ratchet effects, i.e. a rectified transport in absence of a 
net bias force and dissipation. Different ratchet models 
based on optical lattices, which use ultracold atoms as 
ideal systems, have been proposed [14-19], and detailed 
studies have investigated the relation between rectifica- 
tion and broken space-time symmetries [15,17], the ef- 
fects of atomic interactions [19], and the possibility to 
transport quantum information [16,18]. In optics, the 
idea of realizing an 'optical ratchet' has received few at- 
tention to date, and limited to solitary waves in non- 
linear systems [20-22]. In particular, in [20] a setup to 
observe ratchet effects for cavity solitons has been pro- 
posed using coupled waveguide optical resonators with 
an adiabatically-shaken holding pump beam, whereas 
in [21] control of soliton dragging in dynamic optical lat- 
tices has been investigated. 

In this Letter a simple ratchet scheme is proposed in 
a linear optical system corresponding to rectification of 
discretized light refraction in an array of periodically- 
curved zigzag waveguides, which may be realized with 
currently available femtosecond laser writing technology 
[23]. We consider two interleaved arrays of single- mode 
channel waveguides, a primary array and an auxil- 
iary array B^, in the zigzag geometry shown in Fig. 1(a), 
which is similar to the one recently realized in Ref. [23] 
to investigate second-order coupling effects. As opposed 
to Ref. [23] , the propagation constant of modes in waveg- 
uides Bn is assumed to be detuned by a relatively large 
amount cr from the propagation constant of modes in 
waveguides A^ , a condition that can be achieved in prac- 



tice by e.g. lowering the writing velocity of the auxiliary 
array. The optical axis of the waveguides, which lies in 
the (X, Z) plane, is assumed to be periodically-curved 
along the paraxial propagation direction Z with a bend- 
ing profile xq{Z) of period 2zq. In the waveguide ref- 
erence frame z = Z and x = X — xo{Z)^ where the 
arrays appear to be straight, under scalar and parax- 
ial assumptions light propagation at wavelength A is de- 
scribed by the following Schrodinger-like equation (see, 
for instance, [11]) 



ih 



dz 



2n, 



[us - n{x, y)]^p + nsXo{z)xilj (1) 



where h = A/(27r) is the reduced wavelength, Ug is the 
refractive index of the substrate medium, n(x, y) is the 
refractive index profile of the (straight) arrayed struc- 
ture, is the transverse Laplacian, and the dot de- 
notes the derivative with respect to z. As discussed e.g. 
in [11], light propagation in the optical system mimics 
the temporal dynamics of a quantum particle in a bind- 
ing potential rig — n(x, y) driven by an external (inertial) 
force —nsXo{z) with zero mean. In the tight-binding ap- 
proximation, light transfer among the waveguides may 
be described by a set of coupled mode equations for the 
amplitudes and bn of modes trapped in waveguides 
A^ and B^, respectively, which may be cast in the fol- 
lowing form (see, for instance, [24]) 

idn = -naa [(^n+i exp(-2ze) + an-1 exp(2iB)] + (2) 

— i^ah [^n exp(i(7Z — iQ) + hn-i exp(zcrz + zB)] , 
ihn = -tihh [K+i exp(-2z6) + hn-i exp(2iB)] + (3) 

- Kah [dn exp(-zcr2: + zB) + ttn+i exp(-zcr2: - zB)] , 

In Eqs.(2) and (3), Kaa and /^^^ are the coupling rates 
between adjacent waveguides in the primary and auxil- 
iary arrays, respectively, Kah is the cross-coupling rate 
[see Fig. 1(a)], B(z) = 7rnsXo{z)d/ d is the spatial pe- 
riod of the arrays, and a is the propagation constant 
mismatch between waveguides in the two arrays. 
For straight arrays and large mismatch |<j| ^ 
1^ ah •, i^aa-, ^hh^ light dynamics in the two arrays is decou- 
pled and a broad beam, which excites the primary ar- 
ray An at an incidence angle 9 smaller than the Bragg 
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Fig. 1. (Color online) (a) Schematic of a zigzag waveg- 
uide array, composed by a primary array and an aux- 
iliary array B^. (b) Light refraction in a straight waveg- 
uide array, (c) Axis bending profile xq{z) for the obser- 
vation of light refraction rectification. The inset in (c) 
shows a schematic of the periodically-curved array. 
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Fig. 2. (Color online) Rectified refraction (pseudocolor 
map of integrated intensity distribution) in a 65-mm- 
long periodically-curved zigzag waveguide array (up- 
per images), and corresponding behavior expected for a 
straight zigzag array (bottom images), for an input ellip- 
tical beam exciting a few waveguides in the primary ar- 
ray at different tilting angles: (a) ^ = 0, (b) 6> = — 0.56^^, 
and (c) = 0.5^b- The values of other parameters are: 
A = 633 nm, = 1.45, An^ = 0.003, An^ = 0.0032, 
iL' = 3 /im, and (i = 9 /im. 

angle Ob = X/{2d) [see Fig. 1(b)], will diffract and re- 
fract at an angle 6r{0) = 2tvaadsm{7r6/6B) as in an 
ordinary tight-binding linear lattice [4]. Rectification of 
beam refraction, i.e. a net transverse drift of discretized 
beams averaged over all possible incidence angles ^, is 
obviously absent because Or{—0) = —OriO). Even for a 
single array A^, curved to mimic a dc, ac or combined 
dc-ac driving fields, a net transport of wave packets - 
averaged over different initial conditions- is generally un- 



likely, as discussed in Ref. [14]. To achieve a ratchet ef- 
fect, we introduce the auxiliary lattice and assume 
that both arrays are periodically-curved with a bend- 
ing profile ^0(2;) made of a sequence of circular arcs of 
(paraxial) length and constant curvature with alter- 
nating sign xq{z) = ±l/i^, as depicted in Fig. 1(c). Such 
a bending profile corresponds to an ac square- wave force 
—tIsXq^z)^ and to a triangular shape for the function 
Q{z) entering in Eqs.(2) and (3). Note that, since in this 
case the waveguide axis at the entrance plane z = is 
tilted by an angle xo(0) = zq/{2R)^ the incidence angle 
should be now measured with respect to a plane tilted 
by an angle Oc = nsZo/{2R) from the entrance facet [11], 
as shown in the inset of Fig. 1(c). To achieve rectification 
of light refraction in the primary array over one undula- 
tion cycle, the radius of curvature R and spatial period 
zq of alternations are designed to satisfy the conditions 
6(z) = ±a and tZab^o = 7r/2, i.e. 

R = 7rnsd/{X\a\), zq = tt / (2tiah)- (4) 

In the first semi-cycle of undulation, where 6(z) = a, 
neglecting nonresonant terms Eqs.(2) and (3) reduce 
to idn = —i^effbn and ibn = —f^lffCim where A^e// = 
f<iab exp(i(/)) and (/) = a^ — B is a constant phase shift. This 
means that waveguide A^ turns out to be coupled solely 
with waveguide B^, and light transfer between them oc- 
curs like in an ordinary synchronous directional coupler. 
For a propagation length zq equal to one coupling length 
^0 = 'T^/i^Kah)-, light trapped in waveguide A^ is thus 
fully transferred into waveguide B^. At the successive 
propagation semi-cycle, where 6(z) = — cr, Eqs.(2) and 
(3) yield idn+i = -i^effK and ihn = -K^ffan^i with 
^eff = /^a6 exp(z(/)) and ^ = + B, i.e. waveguide B^ 
turns out to be coupled solely with waveguide A^+i. For 
a propagation length zq equal to one coupling length, 
light trapped in waveguide B^ is thus fully transferred 
into waveguide A^+i- Therefore, after a full modulation 
cycle, an arbitrary light distribution in the primary ar- 
ray A^ is shifted, without distortion, by one unit from 
the left to the right. The result is a net drift of refracted 
light, with a locked refraction angle 6r = d/{2zo) inde- 
pendent of the initial incidence angle 6 and beam shape, 
and a suppression of discrete diffraction. Similarly, an 
arbitrary light distribution in the secondary array B^ 
is shifted by one unit but this time from the right to 
the left. As the ratchet effect results from an alternat- 
ing synchronization of couples of waveguides in the array 
{An Bn and Bn ^n+i), different bending profiles 
xq{z), such as sequences of sinuosoidally-curved waveg- 
uides with alternating amplitudes mimicking the driving 
field of Ref. [16], do not generally generate a ratchet 
effect because of incomplete decoupling of the trans- 
port dynamics. As compared to other ratchet schemes in 
tight-binding lattices [14,16,18], our zigzag configuration 
offers the advantage of spatially separating a primary lat- 
tice (in which we want to realize a net transport) from 
an auxiliary one, which is transiently excited. 
The above analytical results are strictly valid within 
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a tight-binding model of Eq.(l) and in the averaging 
Hmit, which neglects weak light transfer among waveg- 
uides due to non-resonant coupling terms. We checked 
the possibility of observing rectification of light refrac- 
tion beyond such approximations by a beam propa- 
gation analysis of Eq.(l). In the simulations, circular 
channel waveguides with a super- Gaussian index profile 
An{x^y) = AnA,B exp{ — [(x^ -\- y'^)/w'^]^}, with radius 
w and index changes An a and Aub in the primary and 
auxiliary arrays, are assumed. The vertical distance dy 
between the two arrays is chosen to be dy = \/3(i/2, 
i.e. the zigzag is defined by an equilateral triangle. As 
an input beam, we typically assumed an elliptical Gaus- 
sian beam, with a tilting angle (in the x direction) 
smaller than the Bragg angle 6b to excite the lowest 
band of the array. The vertical beam position and size 
are adjusted to couple a few waveguides of the primary 
array. Examples of rectified light refraction, observed in 
the curved arrays for a few values of beam incidence an- 
gle ^, are shown in Fig. 2. In the figure, the integrated 
beam intensity distribution J \ip{x^y)\'^dy versus propa- 
gation distance z is shown and compared to the behav- 
ior that one would observe if the arrays were straight, 
i.e. with usual discrete light refraction patterns. For the 
chosen array parameters and wavelength, the coupling 
rates and propagation mismatch are estimated to be 
i^aa = i^bb — i^ab — 1-939 cm~^ and a ~ 16.19 cm~^, 
corresponding to a radius of curvature R = 4 cm and a 
semi-cycle period zq = 8.1 mm according to Eq.(4). The 
numerical results clearly show that, for curved arrays, 
light refraction always occurs at the angle Or = d/{2zo)^ 
regardless of the beam incidence angle ^, and discrete 
diffraction is well suppressed as a result of the controlled 
transport mechanism. Of course beam refraction is re- 
versed when the position of the elliptical input beam is 
shifted to excite waveguides of the auxiliary array. As 
an example. Fig. 3 shows rectification of light refraction 
when the array is excited at the input by an ellip- 
tical Gaussian beam at normal incidence. Note that, as 
compared to Fig. 2 (a), the refraction angle of discretized 
light is now reversed. 

In conclusion, rectification of light refraction in a peri- 
odic photonic structure, which provides an example of an 
Hamiltonian ratchet system, has been proposed, together 
with a possible experimental implementation based on 
periodically-curved zigzag waveguide arrays. 
Author E-mail address: longhi@fisi.polimi.it 
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